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Abstract 

We consider a Lorentz non-invariant dispersion relation for the neutrino, which 
would produce unexpected effects with neutrinos of few eV, exactly where the tri- 
tium beta-decay anomaly is found. We use this anomaly to put bounds on the 
violation of Lorentz invariance. We discuss other consequences of this non-invariant 
dispersion relation in neutrino experiments and high-energy cosmic-ray physics. 
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1 Introduction 

Recent research on the determination of neutrino mass by studying the low- 
energy beta decay spectrum of tritium has produced a best fit value for ml 
which is significantly negative [1]. This unphysical value is caused by an 
anomalous excess of electron events at the end of the spectrum, at about 
20 eV below the end point. The origin of this anomaly is not known. It seems 
clear that there is some systematic effect not taken into account, most prob- 
ably of experimental nature. However, in this letter we want to point that an 
apparent excess of electron events near the end of the spectrum is compatible 
with a certain deviation of the relativistic dispersion relation for the neutrino. 
In this way, the tritium experimental results could be used to put bounds on 
the parameters characterizing this violation of Lorentz symmetry. 
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The idea that Lorentz invariance is not an exact symmetry, but an approx- 
imate one which works extremely well in our low-energy experiments, is not 
new. In the third section of this letter we will review the main proposals for 
modifications of the energy-momentum relativistic dispersion relation. There 
we will argue that a plausible dispersion relation for the neutrino, and only 
for this particle, is the following one: 

El = p u 2 + ml + 2\\ Pu \, (1) 



where |p| means the module of the three- momentum p, m v is the neutrino 
mass, and A is some mass scale, to be determined afterwards. We will consider 
A > in order to have a positive contribution to the energy squared. 



2 Kurie plot 



Let us see how the usual Kurie plot of the tritium beta decay is modified by 
the dispersion relation Eq. (1). The phase-space factor is 

dll = p 2 e dp e pld Pl/ 5(Q-E- E v ) , (2) 



where p e is the momentum of the electron, E is the kinetic energy of the 
electron (E e = m e + E) and Q is the energy available to distribute between 
the neutrino energy and the kinetic energy of the electron. Q is given by 

M-M' = Q + m e , (3) 



M and M' being the masses of the initial and final nuclei. We are neglecting 
here the kinetic energy of the nucleus, which is typically of order 10 _1 eV. The 
Kurie plot is proportional to the function K{E) given by 



K{E) 



J S(Q-E-E v )pldp^ 1/2 



(4) 



For the usual relativistic dispersion relations, E 2 = p e 2 + ml and E 2 = p u 2 + 

m 2 ,, one gets 



K{E) = 



(Q-E)y/(Q-E)*-ml 



■lV2 



(5) 



which is a straight line, (Q — E), in the case m u = 0. 
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If one takes the new dispersion relation for the neutrino, Eq. (1), then the 
function K(E) becomes (from now on, we will use the notation m = m v ) 



K{E) = {Q - E) 1 ' 2 {[(Q - Ef + A 



i2 2 

m 



1/4 



A 



[(Q - E) 2 + A 2 - m 



21 1/4 



(6) 



It is easy to check that the end point of this curve is E = Q — m, just as in the 
standard case Eq. (5). On the other hand, for large Q — E values, the curve 
is very well approximated by a straight line. Considering a point E\ so that 
Q — Ei ^> m, A, we introduce the linear approximation 



K El (E)=K(E l ) + 



dK(E) 



dE 



(E - E l ) , 



(7) 



E=Ei 



i.e., it is the tangent to K(E) at the point E = E\. Expanding Ke { {E) in 
powers of X/(Q — E t ), m/{Q — E[), we get 

K El (E) = (Q-E)-X + ^- )+ ^^ ( E-E l) + ... (8) 



so that, instead of having a straight line of slope —1 ending at E — Q, which 
is the standard case with m = A = 0, we obtain a straight line of slope 



1 A 2 -m 2 

+ 4(g-^) 2 + '-- 



(9) 



which ends at the point 



(10) 



On the other hand, one can calculate the exact value of the slope at the end 
point: 



dK(E) 



dE 



E=Q-m 



m \ 3/2 



(11) 



Two cases are clearly distinguished: A > m and A < m. These are shown in 
Fig. 1. We see that near the end of the spectrum, the curve K(E) is above the 
linear approximation Ke { (E) when A > m, which corresponds to an apparent 
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Fig. 1. The function K{E) for the two cases A > m and A < m, together with 
their linear approximation (dotted line). The straight line Q — E, corresponding to 
A = m = is also shown. Units in this figure are eV. 

excess of electrons at high energies. Indeed it is only apparent, because the 
curve lies always below the corresponding curve of a relativistic dispersion 
relation for a massless neutrino, which is also indicated in the figure. In the 
A < m case, we get the oposite situation: the effect due to the neutrino mass 
dominates over the A term (responsible for the "apparent excess" ) and we get 
a reduction on the number of electrons at high energies. 

The tritium anomaly consists in an excess of electron events at high energies, 
where "excess" here means that the data stay above the straight line which is 
the linear approximation to the curve at low energies. A possible measure of 
the anomaly could be given by using the following quantity 



so that Ra > 1 indicates an apparent excess of electrons at high energies. 
This quantity could be measured experimentally, and a comparison with the 
theoretical prediction would put bounds on the values of the parameters A, m. 
We will show now the orders of magnitude of these bounds in a hypothetical 
but realistic example. 



Ra = 4 



fi 



(g°_ A/2 dEK(E) 
fg°_ A dEK(E) 



(12) 
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Fig. 2. Ra as a function of A/2A for three different values of a: starting from the 
lower curve, a = 0.3 (A = 1.195m), a = 0.75 (A = 2m), a = 0.9375 (A = Am) and 
a = 1 (m/A = 0). The dashed line Ra = 1.0015 is a hypothetical experimental 
bound on Ra (see text). 



Introducing the variable x 
rewritten as follows: 



(Q — E)/\, and taking E ~ Q — A, i?A is 



j*! 2A da; •x/x 



d/4 



far + a) 



-1/4 



A A dx [(X 2 + «)V4 _ ( x 2 + a )-l/4] 



(13) 



where a = 1 — m 2 /A 2 . Fig. 2 shows -Ra as a function of A/2A for different 
values of a, that is, for different values of the quotient m/A. 

Experimentally, one could put a bound on the excess of electron events, that 
is, write -Ra < 1 + ca for a certain value of A, which should be chosen in an 
appropiate way. A should be large enough, so that all the anomaly is contained 
in the region E — A/2 < E < E , but it should not be too large, because 
in that case -Ra — > 1. Let us do an estimate of orders of magnitude. From 
Table 1 of Ref. [2], we see that the fit which gives m 2 < 0, signal of the 
anomaly, is stable and has a good x 2 /d.o.f. in a region between 200 and 400 
eV before the end point, which is Eq ~ 18570 eV. Let us take A/2 = 200 
eV. An approximate value of -Ra comes then from Fig. 2 of Ref. [2], which 
contains the experimental data for the Kurie plot. -Ra is approximately given 
by 4(s + t)/(As + t), where s is the area of a triangle of base and height A/2 
(the straight line of the Kurie plot has slope —1), and t is the area delimited by 
the data above the straight line, spread in a range of energies of around 20 eV, 
and the line itself. A rough estimate of this area gives us a conservative value 
of eA = 1-5 • 10~ 3 . For the four curves of Fig. 2, corresponding to four different 
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values of the quotient m/A, we obtain bounds for A/2A, which, recalling that 
we have taken A/2 = 200 eV, produce bounds for A which go from A < 12.5 
eV to A < 6.5 eV for the extreme curves. Again, this translates into bounds 
for the neutrino mass, depending on the value of a. For A = 2m, the bound is 
m < 3.8 eV, and for A = 4m, we get m < 1.7 eV. 

We see that the tritium anomaly puts bounds for (A,m) in a very interesting 
range (a few eV). Of course, a detailed analysis of the experimental data 
would give finer bounds to this kind of Lorentz invariance violation. Even 
more interesting would be a possible (future) experimental bound of the type 
i?A > 1 + £A) with 6a > (that is, the confirmation of the anomaly), which 
would give a lower bound for A, showing the presence of Lorentz invariance 
violation effects in the tritium beta decay. 



3 Lorentz invariance violations 

Special relativity and Lorentz invariance are at the base of our low-energy 
effective theories. However, it may be possible that these are low-energy sym- 
metries of a larger theory that do not need to be Lorentz invariant. Several 
attempts have been made to question Lorentz invariance [3] and put bounds 
on possible violations (see e.g. [4]). 

One way to explore the potentially observable effects of departures from exact 
Lorentz invariance is to consider the consequences in low-energy processes of 
possible extensions of the Lorentz-invariant particle energy-momentum rela- 
tion compatible with translational and rotational invariance in a "preferred 
frame" . A first possibility is given by the Lorentz- violating class of dispersion 
relations [5] 



where M is the natural mass scale of the Lorentz non-invariant fundamental 
theory. An analysis of cosmic ray processes, whose thresholds are drastically 
changed by the new dispersion relation, puts severe bounds on the scale of the 
Lorentz violation: M has to be several orders of magnitude larger than the 
Planck mass scale [6]. 

Another possible extension of the energy-momentum relation, coming from 
the introduction of a rotationally invariant two-derivative term in the free 
Lagrangian, is [4] 




(14) 



E 2 = p 2 + m 2 + ep 2 , 



(15) 
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where e is a small coefficient which fixes the maximal attainable velocity of 
each particle (v 2 = 1 + e). Differences among maximal attainable velocities of 
different particles lead to abrupt effects when the dimensionless ratio ep 2 /m 2 is 
of order unity. The precise tests of special relativity give very strong constraints 
on this type of extension [4] (e < 1CT 23 ). 

A discussion of departures from exact Lorentz invariance in terms of modi- 
fications of the energy-momentum relation leads to consider a third possible 
extension of the dispersion relation with a term linear in \p\ [7], 

E 2 = p 2 + m 2 + 2X\p\ . (16) 

The additional term dominates over the standard kinetic term (p 2 ) when |p| < 
2A and then the nonrelativistic kinematics is drastically changed. Therefore 
this type of generalized dispersion relation has to be excluded, except just for 
one case. The neutrino has two characteristic properties: it has a very small 
mass, and it interacts only weakly. As a result of this combination, we have not 
any experimental result on its nonrelativistic physics. Therefore, the presence 
of Lorentz invariance violations affecting the nonrelativistic limit cannot be 
excluded a priori in the neutrino case. 

One possible way to incorporate a departure of Lorentz invariance affecting 
only to the neutrinos is to assume that the presence of a linear term in the 
dispersion relation is due to a new interaction which acts as a messenger of 
the Lorentz non-invariance at high energies. Once again, this is not a weird 
assumption: the introduction of new interactions is a general practice in the 
different attempts to explain the smallness of neutrino masses [8]. 

In conclusion, a dispersion relation of the form (1) can be considered for the 
neutrino at low energies. We still assume the existence of further Lorentz 
noninvariant terms, as those contained in Eqs. (14) and (15), for the neutrino 
as well as for all other particles. In the rest of the letter, we turn our attention 
to other implications of Eq. (1). 



4 Other implications of the new dispersion relation 

4-1 Neutrino oscillations 

Let us first concentrate on the influence on a characteristic low-energy phe- 
nomenon for the neutrino: flavour oscillations. In the case of a Lorentz invari- 
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ance violation independent of flavour, that is, 

£ 2 = p 2 + m 2 + 2A|p| , (17) 



one gets the result 



2 2 

E i -E j = ^—^ + ... (18) 
Z \P\ 



so that there is not any footprint of Lorentz noninvariance in neutrino oscil- 
lations. If we admit a possible dependence of A on the flavour, we get 

(A|-A 2 )-(m 2 -m 2 ) , x 

El -E 3 = \- \ 3 - ^ 2i + ... (19) 



and the oscillation probability becomes 



P{vi -> Uj) « sin 2 2# sin 2 L(Km) 



1.27)- 



m 2 )(eV 2 ) 



£(GeV) 



+(3.54)10 9 (A i - A,)(eV) 



(20) 



which means an energy independent probability (which goes against the ex- 
perimental observations in solar neutrino oscillations) unless 



| A; - Xj\ < 10 1S ( 



eV) , 



(21) 



where we have used that L & = 10 8 Km, and A 2 — A 2 mf — m 2 . There will 
be no footprint of Lorentz invariance violations in the case of atmospheric 
neutrinos, since L atm ~ 1O _4 L , nor in accelerator experiments, with much 
smaller lengths. The conclusion is that the linear term in \p\ in the generalized 
dispersion relation (1) has to be flavour- independent. 



4-2 Contribution to the energy density of the Universe 



The bounds on neutrino masses from their contribution to the energy density 
of the Universe [9] are not affected by the presence of Lorentz invariance 
violations, since they only depend on the minimum energy of the neutrino, 
which is still E min = m (for p — 0), independently of A. 
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4-3 Neutrinos in astrophysics 



The spread of arrival times of the neutrinos from SN1987A, coupled with the 
measured neutrino energies, provides a simple time-of-flight limit on m Ue . 

From Eq. (1) one has for the neutrino velocity 

f>=^ = , 1 (22) 
' \P\ 2 (i+VIPI) 2 



Since |p| 3> A, the limit [9] (23 eV) can be taken as an upper limit for 
VA 2 — m? which is not far away from the range of parameters suggested by the 
tritium beta-decay anomaly. Then, neutrinos from supernovas can be a good 
place to look for footprints of Lorentz invariance violations. In fact, we see 
another implication from Eq. (22) if A > m: neutrinos travel faster than light, 
so that those carrying more energy will be the latest to arrive, in contrast with 
the case of a relativistic dispersion relation. 



4-4 Consequences in particle and cosmic ray physics 



A striking consequence of the presence of the A-term in the generalized neu- 
trino dispersion relation Eq. (1) is the kinematic prohibition of reactions in- 
volving neutrinos at a certain energy. Let us consider ordinary neutron decay: 
n — > p + e~ + v. At large momentum, the A-term, which is proportional to 
represents a large contribution to the energy balance, so that the process 
might be forbidden for neutrons of sufficiently high momentum. This turns 
out to be the case: it is easy to see that for X/\p\ <C 1 and \p\ ^> m n ,m p , the 
total energy of the final state is bounded from below, 



E p + E e + E p > \p\ + A + 



2jpT 



(23) 



which implies that neutrons with momentum 



P > 



mz 



[m. 



m e 



2A 



'1.47-10 

2A(eV) 



15 " 



(eV) 



(24) 
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are stable particles. A similar conclusion is obtained in the case of pions: the 
desintegration n~ — > /i~ + ^ is forbidden if the pion has a momentum 



P > 



to; 



to; 



2A 



"8.3 • 10 15 " 
_ 2A(eV) _ 



(eV) 



(25) 



This means that pions and neutrons with these energies can form part of 
cosmic rays, since they are stable particles if the neutrino energy-momentum 
relation is given by Eq. (1). Therefore, cosmic ray physics might be drastically 
affected at high energies, of order E > 10 15 -10 16 eV. In particular, the presence 
of stable neutrons and pions in cosmic rays at such energies might contribute 
to avoiding the well-known GZK cutoff [10]. 



5 Conclusions 



It is surprising that a very simple extension of the Lorentz invariant disper- 
sion relation for the neutrino, consistent with all the constraints from neutrino 
physics, can affect phenomena of such different energy ranges. It could be be- 
hind the tritium beta-decay anomaly, and also lead to a drastic change in the 
composition of cosmic rays beyond 10 16 eV. It seems worthwhile to explore 
in more detail all the consequences of this violation of Lorentz invariance and 
possible models of the Lorentz non-invariant physics at high energies incorpo- 
rating the extended neutrino dispersion relation considered in this letter as a 
low-energy remnant. 



Acknowledgements 



We are grateful to J.L. Alonso for collaboration at early stages of this work, 
to D.E. Groom for useful comments on the tritium anomaly, to A. Grillo, R. 
Aloisio and A. Galante for discussions on observable effects of violations of 
Lorentz invariance, and to J.M. Uson. The work of JMC was supported by 
EU TMR program ERBFMRX-CT97-0122 and the work of JLC by CICYT 
(Spain) project AEN-97-1680. 



References 



[1] V.M Lobashev et al., Phys. Lett. B 460 (1999) 227; Ch. Weinheimer et al, Phys 
Lett. B 460 (1999) 219. 



10 



[2] A.I. Belesev et al, Phys. Lett. B 350 (1995) 263. 

[3] D. Colladay and V.A. Kostelecky, Phys Rev. D 55 (1997) 6760; Phys. Rev. D 
58 (1998) 116002 and earlier references therein. 

[4] S. Coleman and S.L. Glashow, Phys. Rev. D 59 (1999) 116008. 

[5] G. Amelino-Camelia, J. Ellis, N.E. Mavromatos, D.V. Nanopoulos, S. Sarkar, 



Nature 393 (1998) 763; G. Amelino-Camelia, kr-qc/9910089 



[6] R. Aloisio, P. Blasi, P.L. Ghia, and A. Grillo, astro-ph/000125£ , to appear in 
Phys. Rev. D. 

[7] In a recent work, J. Alfaro, HA. Morales- Tecotl and L.F. Urrutia, Phys. Rev. 
Lett. 84 (2000) 2318, it has been shown that a linear term in the neutrino 
dispersion relation appears in the framework of loop quantum gravity. 

[8] See for instance, R.N. Mohapatra in " Current Aspects of Neutrino Physics", 
ed. by D. Cadwell (Springer- Verlag) (2000) and references therein. 

[9] Particle Data Group, Eur. Phys. Jour. C3 (1998) 1. 

[10] K. Greisen, Phys. Rev. Lett. 16 (1966) 748; G.T. Zatsepin and V.A. Kuzmin, 
Zh. Eksp. Teor. Fiz. Pis'ma Red. 4 (1966) 414 [JETP Lett. 4 (1966) 78]. 



11 



